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In the last couple of years, there has been big progress in finite temperature QCD on the lattice. 
Large-scale dynamical simulations of 2+1 flavor QCD with various improved staggered quark 
actions have been started to produce results for various thermodynamic quantities which are ex- 
trapolated to the continuum limit at around physical quark masses, and thus are capable for a 
direct comparison with experiment. At the same time, the theoretical uneasiness with staggered- 
type lattice quarks motivated several groups to accelerate studies with Wilson-type quarks and 
lattice chiral quarks. In this review, I discuss these important developments in finite temperature 
QCD made in the past year. 
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1. Introduction 

The transition from the conventional hadronic matter to the quark-gluon plasma (QGP) is one 
of the most drastic predictions of QCD. The rapid progress of heavy ion collision experiments urge 
us towards quantitatively reliable determination of properties of hot quark matter directly from the 
first principles of QCD. In this note, I review the recent status of studies on the nature of the quark 
matter at finite temperatures on the lattice. 

Study of finite temperature QCD often requires simulations on coarse lattices. It is thus impor- 
tant to keep in mind possible caveats from lattice artifacts. The choice of the lattice quark action 
has significant implications for the nature of lattice artifacts. Therefore, before discussing individ- 
ual topics, I first comment in Sect. 2 on the nature of various lattice quarks adopted in the studies 
of finite temperature QCD, as well as possible caveats from them. 

Because the s quark mass is comparable with the QCD scale Aqcd and thus the transition 
temperature Tc, it is important to incorporate the dynamical effects of the s quark for realistic 
predictions about QGP. In Sect. 3, 1 discuss the status of our understanding of the phase structure 
of QCD with 2-1-1 flavors of quarks. There was a conflict about the value of Tc with improved 
staggered quarks. This year, the discrepancy is mostly resolved through intensive simulations on 
fine and large lattices. From these developments, we learn the importance of a good control of 
lattice artifacts in order to achieve a consistent prediction for experiment. Scaling is a powerful 
tool to discriminate the nature of the transition. I discuss recent progress in the studies on the 
scaling properties of the QCD transition. 

Section 4 is devoted to the status of the equation of state (EOS) in 2-1-1 flavor QCD. The scope 
of finite temperature QCD is quite wide. Some of important developments in various directions 
are discussed in Sect. 5. Studies of finite density QCD will be summarized by S. Gupta in these 
proceedings [1]. For recent reviews, see [2, 3, 4]. 

2. Quarks on the lattice 

2.1 Staggered-type quarks 

Most quantitative outputs of finite temperature QCD have been produced with staggered-type 
lattice quarks. We have good reasons to adopt these quarks. Namely (i) they preserve a kind of 
chiral symmetry (U(l) symmetry under a taste-chiral transformation) on the lattice, which protects 
the location of the chiral limit against quantum fluctuations, and (ii) they are relatively cheap to 
simulate. The latter is important for studies of finite temperature QCD which usually require a 
systematic survey in a wide range of the parameter space. 

On the other hand, the staggered quarks lead to four copies ("tastes") of fermions for each 
flavor in the continuum limit. To remove unwanted three tastes, we usually adopt the "fourth 
root procedure", i.e. we replace the quark determinant by its fourth root in the path-integration 
procedures. Since this makes the theory non-local and non-unitary, many theoretical issues emerge 
[5, 6, 7, 8]. In particular, the universality arguments become fragile due to non-local interactions, 
and thus, e.g. we are not fully sure if the staggered quarks have the continuum limit with desired 
number of flavors. So fai^, no apparent obstacles about the validity of rooted staggered quai^ks seem 
to exist, provided that the continuum extrapolation is done prior to the chiral extrapolation [5, 8]. 
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It is out of the scope of this review to go into these discussions. Therefore, in the followings, I 
simply assume that the continuum extrapolation is (will be) done first and that the staggered quarks 
do have the continuum limit which is in the desired universality class of the QCD with arbitrary 
number of flavors. 

Even with these assumptions, a couple of issues remain. One of them is the problem of taste 
violation. The symmetry among tastes in the continuum limit is explicitly violated at finite lattice 
spacings. This introduces a systematic uncertainty in the identification of flavors. Due to the taste 
violation, for example, we observe different masses for the same hadron depending on the combi- 
nation of tastes for the hadron operator. The largest taste violation is known to appear in the pions 
[9]. With the U(l) taste-chiral symmetry preserved on the lattice, one combination of tastes, the 
pseudo Nambu-Goldstone (pNG) pion, leads to the vanishing mass in the chiral limit, but other 
15 combinations remain heavy. Although the lightest pNG pion is customary treated as the pion, 
heavier pions do contribute to loop diagrams and thus induce lattice artifacts. To reduce the mag- 
nitude of lattice artifacts, various improved actions have been proposed, such as asqtad [10], HYP 
[11], p4 [12], stout [13], HISQ [14], etc. Each of them may be originally designed to improve 
different quantities, such as the free quark dispersion relation, perturbative propagators etc., but re- 
duces the taste violation too directly or indirectly. Note that, unlike the case of hadron propagators, 
loop diagrams including all tastes are the leading contributions in many thermodynamic quantities. 
Therefore, a good control of the taste violation is important in finite temperature QCD. 

Another problem which is worrisome in particular in finite temperature QCD is if we can 
expect scaling properties at finite lattice spacings in spite of non-local interactions. This matters, 
e.g. when we study scaling properties around a second order chiral transition or near a critical end 
point of a first order ti'ansition line. I know no theoretical arguments supporting applicability of 
universality arguments with rooted staggered quarks on finite lattices. 

I come back to these issues later. 

2.2 Wilson-type quarks 

A more conservative approach is to adopt Wilson-type lattice quarks with which the existence 
of the continuum limit is guaranteed for any number of flavors, at the price of explicit violation of 
the chiral symmetry at finite lattice spacings and more computational costs near the massless limit. 
In the studies with Wilson-type quarks, we need to keep in mind that the phase structure is not 
simple at finite lattice spacings due to the explicit violation of the chiral symmetry. In 1984, Sinya 
Aoki predicted the existence of a second order phase transition line in the phase diagram to explain 
massless pions without resorting to chiral symmetry [15]. Beyond the second order transition line, 
we have the "Aoki phase" in which a parity-flavor symmetry is spontaneously broken. Fourteen 
years later, Sharpe and Singleton pointed out that the chiral violation of Wilson quarks induces 
an additional term in the chiral effective theory, and the phase structure at finite lattice spacings 
changes depending on the sign of the coiTcsponding effective coupling C2 [16]. Namely, when C2 
is positive, we have the second order transition as predicted by Aoki. But, when C2 is negative, we 
instead have a first order transition. 

The phase structure becomes more transparent when we extend the action to include a mass 
twisting term [17], S^™ = a''^\\f(x) [Dw + nio + 'Ys^sMo] Wi^)'' where Dw is the Wilson quai^k ker- 
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Figure 1: (Left, Center) Phase structure for twisted Wilson-type quarks at finite lattice spacing as a function 
of the untwisted bare mass parameter mo and the twisting parameter /io, based on a prediction by Ref. [18]. 
(Right) Suggested phase structure for twisted Wilson-type quarks at finite Nt [21]. The vertical axis is the 
hopping parameter. Possible first order cuts in the jU direction are hidden in the regions marked by "?". 



nel and T3 acts in the isospin space. Sharpe and Wu [18] studied the chiral effective theory for 
twisted quarks and found a phase structure schematically plotted in the left and center panels of 
Fig. 1. The untwisted theories {jiQ = 0) are on the horizontal axis. For the case of C2 > 0, we have 
the first order cut corresponding to the Aoki phase on the horizontal axis. When we decrease the 
mass mo along the horizontal axis, we hit the second order boundary of the Aoki phase, as shown 
by a red arrow in the figure. When C2 < 0, on the other hand, the cut appears along the vertical axis 
and we encounter a first order transition when we decrease iiiq. 

Although the sign of C2 should depend on the details of the lattice action, various gauge actions 
combined with twisted standard Wilson quark action are known to have the first order transition at 
large j8 . In this case, we have to avoid the first order region to approach the continuum Umit. While 
the width of the first order region is expected to shrink rapidly towards the continuum hmit, a region 
with sufficiently light quarks may exist only at small lattice spacings. On the other hand, the PACS- 
CS Collaboration have studied QCD with non-perturbatively improved clover quarks coupled to the 
RG-improved Iwasaki glue, and succeeded in decreasing the light quark mass down to the physical 
point without encountering a first order transition [19]. Therefore, for this combination of the 
actions, either C2 > is realized or the first order transition region is located at quark masses lighter 
than the physical light quark mass. The phase structure with twisted improved quarks is not well 
clarified yet [20]. 

At finite temperature, Creutz [21] studied a chiral effective model for twisted Wilson quarks 
and suggested that the deconfinement transition surface in the extended parameter space forms a 
deformed cone-like shape, as shown in the right panel of Fig. 1 . The tmfT Collaboration studied 
twisted standard Wilson quarks, and obtained signals consistent with the predicted cone-shaped 
deconfinement transition surface [22, 23]. Further studies are needed to clarify the full phase 
structure for twisted Wilson-type quarks. 

2.3 Lattice chiral quarks 

The most attractive but currently the most expensive way to push forward is to adopt a chiral 
lattice quark action. Domain-wall (DW) quarks [24] are formulated on a 5-dimensional space-time 
and designed such that the chiral symmetry is recovered at finite lattice spacings in the limit L, — )• 00 
where Lg is the lattice size in the 5th dimension. In actual simulations, however, we have to keep 
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finite. Tliis causes a residual cliiral violation, wliicli is conventionally measured by the magnitude 
of the additive renormalization term m^^s to quark masses near the continuum limit, mjes is called 
the residual quark mass and is known to decrease exponentially with for small L,, but then 
levels off into a much slower decrease proportional to 1 /Lg due to the mobility edge phenomenon 
[25]. The requirement of large computational resources forces us to work on lattices coarser than 
other lattice quarks and with small L^. Because DW quarks require smooth gauge configurations, 
coarse lattices may introduce serious chiral violations. A study of finite temperature QCD with 
DW quarks was started about 10 years ago for A^^ = 2 on a coarse lattice (8^ x 4) with L, 8 [26]. 
Although the chiral violations seem to be controled by mres, the values of /Wres were too large to 
discuss physical properties. 

This year, a new study of 2+1 flavor QCD 
adopting the RG-improved Iwasaki gauge action 
and much larger « 32 on a finer lattice (16^ x 8) 
has been published [27]. Both the improvement 
of gauge action and the adoption of a larger are 
effective to reduce mres- Accordingly, their results 
for chiral susceptibilities around the thermal tran- 
sition point show an effective quark mass depen- 
dence qualitatively consistent with expectations. 
However, even with these improvements, the re- 
sulting wires turned out to be not small enough: 
Wres ~ 0.008 is larger than the bare light quark 
mass mi = 0.003 needed to achieve the effective 
light to strange quark mass ratio (wi/ +mies) / (wJ.v + 
wires) ~ 0.025. The HotQCD Collaboration is now 
extending the study adopting an even larger Lg of 
96, and also with an improved action [28] dedicated for DW quarks. See Fig. 2 [29]. 

Finally, Cossu reported at the conference the status of a finite temperature study that JLQCD 
collaboration has started using an overlap quark at a fixed topology [31]. 

3. Phase structure of QCD at finite temperatures 

Figure 3 summarizes the current wisdom about the order of the finite temperature transition 
in 2+ 1 flavor QCD as a function of the quark masses, based on the studies of lattice simulations 
and effective models. When the u, d and s quarks are all sufficiently light or sufficiently heavy, we 
expect the transition to be of first order. At intermediate values of the quark masses, the "transition" 
will be actually an analytic crossover. On the boundaries of the first order regions, we expect second 
order transitions. Associated critical scaUng around there will have characteristics universal to the 
con^esponding effective models. 

Simulations using improved staggered quarks suggest that the physical point, identified by the 
pNG pion mass and other hadron masses, is located in the crossover region (see e.g. [32]). Now, 
the main challenges are to compute the detailed properties of the transition at the physical point 
with a high precision. The transition temperature Tc is one of the most important quantities for 




T (MeV) 



Figure 2: Disconnected chiral susceptibility ob- 
tained with 2+1 flavors of DW quarks (red square) 
[27], compared with the results of asqtad and 
HISQ quarks [30] Preliminary results for = 96 
(black circles) and with an improved gauge action 
(blue downward triangles) [29] are also shown. 
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Figure 3: Order of the finite temperature transition in 2+1 flavor QCD as a function of the degenerate u 
and d quark mass m^d and the s quark mass mj. (Left) The standard scenario with the second order chiral 
transition for two-flavor QCD. (Right) An alternative scenario when the two-flavor chiral transition is first 
order. 

experimental investigations of QGP. Estimation of Tc in 2-1-1 flavor QCD has been made based on 
large-scale simulations using various improved staggered quarks. However, there has been a big 
discrepancy in the values of among different groups for more than five years. This year, the main 
part of the discrepancy has been removed. 

The nature of the transition in the chiral limit of two-flavor QCD (the upper left edge of the 
figure) has significant implications for the nature of the transition at the physical point too. The 
left panel of Fig. 3 summarizes the standard scenario in which the chiral transition of two-flavor 
QCD is second order in the universality class of the 0(4) Heisenberg model [33]. In this case, 
because the chiral transition of thi^ee-flavor QCD is of first order, we have a tricritical point on 
the left edge of the figure (m,„/ = 0) where the order of the transition changes from the second 
order to the first order. Depending on the location of the tricritical point relative to the physical 
point, the universality class dominating the parameter dependence around the physical point will 
be different. The right panel of Fig. 3 shows an alternative scenario in which the chiral transition of 
two-flavor QCD is first order. In this case, we have no tricritical point and thus no regions for the 
0(4) universality class. A distinction between the two scenarios is important for studies at finite 
densities too. Although the majority view the standard scenario as more probable, the nature of the 
two-flavor chiral transition was not fully fixed. This year, we had some advances. 

In this section, I discuss these developments. 

3.1 Transition temperature 

In 2005, the MILC Collaboration obtained Tc = 169(12) (4) MeV in the combined chiral and 
continuum limit from a measurement of a chiral susceptibility in 2-1-1 flavor QCD with asqtad 
quarks and the one-loop Symanzik glues on A', = 4-8 lattices [34], where the scale was set by ri 
and the 0(4) critical exponent was adopted in the chiral extrapolation. In 2006, the Wuppertal- 
Budapest Collaboration has published their values based on a study of the 2-1-1 flavor QCD with a 
stout-link improved staggered quarks coupled to the tree-level Symanzik glues [35]. Carrying out 
a chiral extrapolation to the physical point and a continuum extrapolation using A', = 6-10 lattices. 
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they found Tc = 151(3)(3) MeV from the susceptibility of the chiral condensate, = 175(2)(4) 
MeV from the s quark number susceptibility, and Tc = 176(3) (4) MeV from the inflection point 
of the Polyakov-loop, where the scale was set by fx- In the same year, Cheng et al. (let me call 
them as "BNL-Bielefeld Collaboration") pubhshed T^ = 192(7) (4) MeV at the physical point in 
the continuum limit, based on their study on A/, = 4 and 6 lattices using the p4fat3 staggered quark 
action and tree-level Symanzik gauge action [36]. Their Tc is an average of estimates from the 
Polyakov-loop susceptibiUty and the chiral susceptibility, and the systematic error includes the 
difference between them. Later in 2008, the HotQCD Collaboration (which is a merger of the 
BNL-Bielefeld Collaboration and the MILC Collaboration) compared the results of p4 and asqtad 
actions, and found that the two staggered quarks lead to roughly consistent results for the location 
of the transition on finite lattices [37, 38]. 
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Figure 4: Chiral observables in 2+1 flavor QCD with improved staggered quarks. The HISQ, p4 and asqtad 
data are obtained at the bare quark mass ratio niud /nis= 0.05 with around the physical value, using a scale 
fixed by ro . The physical point is estimated to be niud/ms ~ 0.037. The stout data, which is obtained at the 
physical point using a fK scale [40], is shifted in T to correct the difference in the scale setting. Note that this 
procedure causes a slight deviation from the physical point at finite lattice spacings. (Left) Subtracted chiral 
condensate (difference of scaled light quark and s quark condensates to remove a divergent renormalization 
factor) from the HISQ, p4, asqtad and stout actions [30]. (Center) Comparison of HISQ and asqtad results 
for the disconnected part of the chiral susceptibility [30]. (Right) Comparison of HISQ and stout results for 
the renormalized chiral susceptibility in which the T = Q contribution is subtracted [42]. 



Because we have an analytic crossover instead of a singular phase transition, the value of Tc 
may depend on the choice of observable to define it. Details of the analyses, such as the scale 
setting and the definition of the line of constant physics, are different among the two groups too. 
However, it was shown that the discrepancy exists even when we adopt the same observable and 
the same scale convention [39]. See the left panel of Fig. 4. We find that the whole temperature 
dependence is shifted depending on the choice of the action. 

Both groups have pushed forward their studies to clarify the origin of the discrepancy. The 
Wuppertal-Budapest Collaboration extended the study to finer lattices (A^, up to 16) and concluded 
that their values remain essentially unchanged [39, 40]. Their final values are Tc = 147(2)(3) 
MeV from the susceptibility of the chiral condensate and Tc = 165 (2) (4) MeV from the s quark 
number susceptibility [40]. On the other hand, the HotQCD Collaboration noted that the transition 
temperature from p4 and asqtad quarks show rather rapid shift of —5 to —7 MeV on finite lattices 
when the lattice spacing is decreased from A', = 6 to 8 [38], and another shift of about —5 MeV 
when the bare light quark mass is decreased from m„d/ms = 0.1 to 0.05, that corresponds to a 



7 



Finite Temperature QCD on the Lattice - Status 2010 



Kazuyuki Kanaya 



decrease of the pNG pion mass from about 220 MeV to 160 MeV [41]. At the conference, they 
presented a preliminary result = 164(6) MeV at the physical point in the continuum limit, from 
their study of disconnected chiral susceptibility using asqtad quarks on lattices up to A', = 12 [30]. 
At the same time, the Hot QCD Collaboration has presented first results of their project on QCD 
thermodynamics adopting a new improved staggered quark action, the "highly improved staggered 
quark (HISQ)" action [30, 42]. A comparison of these actions is shown in Fig. 4. 

Consulting the recent values of Tc and the location of the peaks in the chiral susceptibilities 
shown in Fig. 4, we find that the discrepancy in Tc is now mostly removed: All the actions converge 
towards T^- ^ 145-165 MeV. 

What was the origin of the discrepancy? The main difference among the improved staggered 
quark actions seems to be the magnitude of the remaining taste violation: When we adjust the bare 
quark mass to have ^135 MeV at lattice spacings corresponding toT ^ 170 MeV on A', = 8 
lattices, the masses for the heavy pions are around 400-600 MeV (asqtad), 300-500 MeV (stout) 
and 200^00 MeV (HISQ), respectively. The p4 is slightly worse than the asqtad. Therefore, the 
p4 and asqtad quarks suffer from severer lattice artifacts than the stout and HISQ quarks. Note that 
the simulations with stout quarks were performed on finer lattices. On an A', = 12 lattice, the heavy 
pion masses with the stout quark action will be around 200-350 MeV at T 170 MeV. 

The overall effects of light and heavy pions in loop 
diagrams may be measured by a root mean squared 
(RMS) mass of pions m^^^ [3]. In Fig. 5, the Wuppertal- 
Budapest Collaboration demonstrated that the asqtad 
data can be reproduced by the stout quark action by 
adjusting [40]. The open rectangles are results 

of stout quark action at rn^^ = 414 MeV [40], and 
the red curve is the result of asqtad quark action at 
m^'^ = 220 MeV [38]. The stout quark mass for open 
rectangle is adjusted to reproduce = 587 MeV of 
the asqtad data atT ^ 135 MeV on A', = 8. According 
to the larger taste violation with the asqtad action, the 
asqtad data at a given m^^^ corresponds to the stout 
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Figure 5: Subtracted chiral condensate with 
stout and asqtad quarks at M = 8 [40] . 



data at a larger . This explains the larger lattice artifacts in the p4 and asqtad results at A^, < 



While the major part of the discrepancy is explained by the taste violation, small disagreements 
remain. In particular, the sensitivity of the determination of T^ on the operator, observed with the 
stout quark, is not confirmed by other quarks yet. 

In all of these studies, the physical point is identified by the pNG pion mass. The success 
of the RMS pion mass, however, suggests that it is more appropriate to consult the RMS hadron 
masses to judge the location of the physical point for thermodynamic quantities. This will make 
the physical point closer to the chiral limit and the values of T^ even smaller on finite lattices. 

It is highly desirable to confirm the results with other types of lattice quarks. Unfortunately, 
the current status with Wilson-type quarks is quantitatively much behind. Recent updates of T^ in 
two-flavor QCD ai^e given in Refs. [23, 43, 44]. The first rough estimate of T^ with 2-1-1 flavors of 
DW quarks is given in [27]. 
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3.2 Chiral scaling 

At r = 0, QCD in the chiral limit has the global chiral-flavor symmetry \]{Nf)r x U(A^f )l 
naively, where is the number of flavors. Among them, U(1)a is explicitly violated by the 
anomaly, and the remaining symmetry is dynamically broken down as SU(A^/r )r x SU(A'^f)l x 
U(l)v — 5- SU(A'^f)v X U(l)v- At sufficiently high temperatures, the broken symmetry will be 
restored. From a study of effective Ginsburg -Landau theory respecting these symmetries [33], the 
chiral restoration transition is predicted to be of first order for Np > 3. For the case of Nf = 2, 
the transition is predicted to be either second order in the universality class of 0(4), or first order 
if the anomaly is weak around the transition point. Clarification of the Nf = 2 chiral transition 
is important for understanding the nature of the transition at the physical point, as summarized in 
Fig. 3. Here, scaling analysis is a powerful tool to discriminate the order of phase transitions. 

3.2.1 0(4) scaling with Wilson-type quarks 




Figure 6: 0(4) scaling fits with in two-flavor QCD with Wilson-type quarks. (Left) Chiral condensate with 
the unimproved Wilson quark action and the Iwasaki gauge action [46]. The chiral condensate is rescaled 
using 0(4) critical exponents, and the curve is the scaling function obtained in an 0(4) spin model. (Center) 
Rescaled chiral condensate with a MF-improved clover quark action and the Iwasaki gauge action [47]. The 
dashed curve is the scaling function obtained in an 0(4) spin model. (Right) 0(4) fit to the pseudocritical 
temperature T^. with a NP-improved clover quark action and the unimproved plaquette gauge action [43]. 

With Wilson-type quarks, the 0(4) scaling expected for the case of a second order chiral transi- 
tion has been observed with various combinations of Wilson-type quark actions and gauge actions, 
since the early years of the investigations. Although the chiral symmetry is explicitly violated on 
the lattice due to the Wilson term, when the additive renormalizations are coiTcctly performed such 
that the chiral violations are subtracted out and the chiral Ward identities are satisfied in the con- 
tinuum limit, the resulting subtracted observables will show the conect chiral behaviors near the 
continuum limit [45] . 

Adopting a chiral condensate defined by the axial vector Ward identity [46], the QCD data are 
shown to reproduce the scaling function of the 0(4) Heisenberg model, with both unimproved and 
clover-improved Wilson quark actions [46, 47, 48]. See the left and center panels of Fig. 6. As 
shown in the right panel of the same figure, a recent study on finer lattices with slightly lighter quark 
masses supports the 0(4) behavior too [43]. These results, together with several observations made 
in these papers, suggest that the transition is continuous in the chiral limit. However, the quai^k 
masses studied are quite heavy yet (niji > 400 MeV). Confirmation with lighter quai^ks on a fine 
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lattice is indispensable. At the conference, Brandt reported the status of their study towards this 
direction [49]. 

3.2.2 Chiral scaling with staggered-type quarks 

Scaling studies with staggered-type quarks have caveats. A^f / 4 is realized through the fourth 
root procedure discussed in Sect. 2.1. Therefore, the symmetry of the system at finite lattice spac- 
ings is that of the original four-taste theory. In the chiral Umit on finite lattices, we thus have the 
0(2) [= U(l)] symmetry for any number of flavors. 

This means that, when the chiral transition is of second order for Np = 2, we should expect 
scaling properties in the universality class of 0(2) spin models, instead of the 0(4). Here, we have 
assumed that the non-locality of rooted staggered quark actions does not invalidate the universality 
arguments at finite lattice spacings too. In any case, the 0(4) universality may appear only when 
we take the continuum limit prior to the chiral extrapolation. We may, however, expect that when 
the chiral transition is continuous on finite lattices, it remains so in the continuum limit. In this 
sense, examination of 0(2) is important, with the caveats mentioned above. 

Numerically, it is not very easy to distinguish between 0(4) and 0(2) scalings with the current 
accuracy of lattice simulations because the critical exponents and the scaling functions in both 
theories are almost identical. 
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Figure 7: (Left) 0(4) scaling test of the chiral condensate in two-flavor QCD with unimproved staggered 
quark action and unimproved plaquette gauge action on A^, = 4-12 lattices. The curves are the 0(4) scaling 
function [51]. (Right) Histogram of chiral condensate in two-flavor QCD with unimproved staggered quark 
action and unimproved plaquette gauge action on an A, = 4 lattice [54]. 

Besides the caveats, the results for the chiral transition with two flavors of staggered quarks 
have been quite confusing. Although the early studies suggest that the transition is continuous 
in the chiral limit, the scaling behaviors observed there show neither 0(2) nor 0(4) [50, 51, 52]. 
The left panel of Fig. 7 is a comparison of QCD data with the 0(4) [w 0(2)] scaling function by 
the MILC Collaboration [51]. 0(2) or 0(4) scaling may appear only on very fine lattices and/or 
at smaller quark masses. (With a modified model including a four-fermi interaction, appearance 
of the 0(2) scaling is reported at vanishing quark masses [53].) In contradiction to these studies, 
the Genova-Pisa group reported indications of a first order transition on an A'j = 4 lattice [54]. A 
possible weakly two-state signal observed by them is shown in the right panel of Fig. 7. 
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Figure 8: Rescaled chiral condensate Mq = ms{\lfY}ud/T'^ in 2+1 flavor QCD with p4 staggered quark 
action and tree-level Symanzik gauge action on A', = 4 [55] (Left) and 8 [56] (Right) lattices. The QCD 
data are fitted to the 0(2) scaling function with an ansatz for deviations from the scaling. 



All of these studies are made using the unimproved staggered quark action and the unimproved 
plaquette gauge action. As discussed in Sec. 3.1, large taste-violation in these studies will make 
the quark masses effectively much heavier than those naively estimated by the pNG pion mass. 

This year, a new study with improved staggered quark action was presented: Adopting the p4 
quark action and the tree-level improved Symanzik gauge action, the BNL-Bielefeld Collaboration 
performed 0(2) and 0(4) scaling fits in 2-1-1 flavor QCD on an A'f = 4 lattice for the bare quark 
mass ratio in the range ni,„//ms = 0.05 down to 0.0125 with around the physical point [55]. 
The quark mass range corresponds to m^'^ s» 150 MeV down to 75 MeV. At the conference, 
preliminary results obtained on an Nt = S lattice were presented too [56]. The absence of spatial 
volume dependence suggests that the simulation points are in the crossover region. For the first 
time with staggered-type quai^ks, 0(2) [or 0(4)] scaling has been observed. At the same time, they 
found substantial deviations from the 0(2) scaling function for m„d/ms > 0.05. In Fig. 8, light 
quark chiral condensate is fitted to the 0(2) scaling function assuming 0(2) critical exponents, 
together with correction terms modeling the deviations from the scaling [55]. Different curves 
represent deviated scaling function for each quark mass. They have tested the 0(4) scaling too. 
However, discrimination of 0(2) from 0(4) was not possible numerically. 

Together with the adoption of improved actions, the substantial deviations from the 0(2) scal- 
ing for niud/ms > 0.05 explain the failures in previous studies. A big impact of this study is the 
0(2) dominance around the physical point. First, this suggests that the chiral transition in two- 
flavor QCD is second order. Second, the s quark mass for the tricritical point may be smaller than 
the physical value. In Fig. 3, 1 have therefore put the physical point slightly higher than the previ- 
ous plots. The second observation is consistent with the fact that the region of first order transition 
around the three-flavor chiral limit has shrunk largely through the introduction of improved actions 
and the adoption of larger A'f in the studies with staggered-type quarks. When the parameter- 
dependence around the physical point is confirmed to be dominated by the 0(4) exponents and the 
0(4) scaling function in the continuum limit, the studies of QGP will be largely accelerated. 

3.3 From the heavy quark limit 

Before closing this section, let us turn our attention to the upper right corner of Fig. 3. In the 
heavy quark limit, the deconfinement transition of SU(3) pure gauge theory is described by the first 
order transition of effective Z(3) Potts model [57]. The phase of Polyakov loop serves the role of 
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Figure 9: (Left) Slope of an effective potential for the plaquette as a function of the hopping paremeter 
(f, obtained in two-flavor QCD with the standard Wilson quark action and the plaquette gauge action on an 
Nt — 4- lattice [62]. (Center) The phase diagram for 2+1 flavor QCD around the heavy quark limit [62]. 
(Right) Trace anomaly in quenched QCD [64] . 



Z(3) spin. Dynamical heavy quarks affect as an external magnetic field to the effective spins and 
break the Z(3) symmetry explicitly. When we decrease the quark mass, the external magnetic field 
becomes stronger and eventually wipes away the first order transition into a crossover. 

In the heavy quark limit, the Z(3) symmetry is spontaneously broken in the high temperature 
phase. The deconfinement transition may thus be described by the percolation of Z(3) clusters. Ex- 
tending a previous study in the pure gauge theory [58], Danzer et al. studied the effect of dynamical 
quarks on the Z(3) clusters defined by the phase of local Polyakov loops [59]. The percolating na- 
ture at high temperatures turned out to be similar to the case of the pure gauge theory. However, 
the temperature dependence is much smoother than the pure gauge theory, in accordance with the 
crossover expected at the physical point. Langelage, Lottini, IVIiura, Nakano and Ohnishi discussed 
about Polyakov loop effective theories in the context of a strong coupling expansion [60, 61]. 

H. Saito reported a study on the fate of the deconfinement transition when the quark mass is 
decreased from infinity, adopting a method using probability distribution of observables combined 
with a reweighting technique [62]. They succeeded in calculating the effective potential, defined 
as the logarithm of the probability distribution for the plaquette, for a wide range of the plaquette 
expectation value in the heavy quark mass limit. They then incorporated the effects of dynamical 
quarks by reweighting combined with the hopping parameter expansion. The left panel of Fig. 9 
is their result for the slope of the effective potential as a function of the hopping parameter K. We 
see that the "S" shape typical for a first order transition gets weaker when K is increased {nicj is 
decreased). Determining the location of the critical point terminating the first order transition, they 
obtained the phase diagram shown in the center panel of Fig. 9. As naively expected, the endpoint 
monotonically shifts to lighter quark mass as we increase the number of flavors. The location of 
the end point is consistent with a previous estimation using an effective Z(3) spin model [63]. 

Borsanyi studied the EOS in quenched QCD up to extremely high temperatures to examine 
the onset of perturbative behavior [64]. Adopting an iterative technique and various improvements, 
they calculated the EOS up to 10000 as shown in the right panel of Fig. 9. They found that their 
result is consistent with a NNLO hard thermal loop perturbation theory above 8 T^. 
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4. Equation of state 

Calculation of the EOS in full QCD is demanding. Besides the basic information such as the 
lattice scale, values of the beta functions, and precise shape and location of the line of constant 
physics (LCP), we need to know the expectation values at T = to renormalize finite temperature 
observables at each simulation point. Conventionally, we vary the coupling parameters along a LCP 
on a lattice with a fixed temporal lattice size A^, to vary the temperature (fixed Nt approach). To 
calculate the EOS non-perturbatively, the integration method is adopted for which precise data at 
lower temperatures are needed [65]. These require a systematic scan in a wide range of parameters 
on both zero and finite temperature lattices. Furthermore, for a continuum extrapolation, we have 
to repeat the calculation at several values of A^,. 

Due to the big demand on the computer power, calculation of the EOS in realistic 2+ 1 fla- 
vor QCD has been made almost exclusively with staggered-type quarks. This year, high preci- 
sion results of the EOS for (almost) physical quark masses, identified using the pNG pion mass, 
were presented adopting different improved staggered quark actions, though lattices are slightly 
coarser than those discussed in Sec. 3.1. The first calculation of the EOS in 2-1-1 flavor QCD with 
Wilson-type quarks has been reported too, in which a new method was developed to beat the higher 
computational demand with Wilson-type quarks. 

4.1 EOS with 2+1 flavors of improved staggered quarks 




140 160 180 200 220 240 260 140 160 180 200 220 240 260 150 200 250 300 350 400 450 500 



Figure 10: (Left, Center) EOS in 2+1 flavor QCD at w 154 MeV (ni,,d/ms = 0.05) obtained on an 
Nt = 8 lattice with p4 staggered quarks [41]. Results at m'^'^ = 220 MeV (mnj/nis = 0.1) are compared. 
(Right) Comparison of the trace anomaly in 2h-1 flavor QCD with HSQ, asqtad and p4 quarks at mud/nis — 
0.05 on A/, = 8 lattices [30]. 

The BNL-Bielefeld Collaboration published the EOS "with physical quark masses", based on 
their simulations using p4 quarks on an A^, = 8 lattice at the bare quark mass ratio m„^//m.v = 0.05 
with the s quark mass fixed near the physical point [41]. Using a scale set by the Sommer scale 
ro = 0.469 fm, their simulation points correspond to m^'^ = 154 MeV, m^^'^ = 486 MeV and 
= 663 MeV with a precision of 3%. This study is an extension of their previous papers on 
the EOS with "almost physical quark masses' ' {mt,d/m, = 0.1 where m^^'^ = 220 MeV) [37, 38]. 
In the left panel of Fig. 10, the trace anomaly ar ound the transition temperature is shown. The shift 
of about —5 MeV from mud/ms = 0.1 mentioned in Sec. 3.1 is visible. Adopting the integration 
method, the EOS shown in the center panel is obtained. Note that the physical point was identified 
by the pNG pion mass, and the chiral and continuum extrapolations are not done yet. 
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Figure 11: (Left) Trace anomaly in 2+1 flavor QCD with stout, asqtad and p4 quarks on lattices Nt>S [66]. 
The stout values are corrected by a tree-level improvement factor. (Center, Right) EOS in 2+1 flavor QCD 
with stout quarks [66]. 



At the conference, Bazavov and Soldner presented preliminary results of the HotQCD Collab- 
oration for the EOS obtained with HISQ and asqtad quarks at m„f//nj4 = 0.05 on Nt = S lattices 
[30]. In the right panel of Fig. 10, the results of the trace anomaly are compared at m„^/m.v = 0.05 
for HSQ, asqtad and p4 quarks. The location of the transition is well consistent with each other, as 
discussed in Sec. 3.1. On the other hand, the values of pressure and energy density are still moving. 
We note that the peak of the trace anomaly around T ~ 200 MeV becomes lower with the adoption 
of a more improved action. 

The Wuppertal-Budapest Collaboration presented the results for the EOS with stout quarks at 
the physical point, obtained on Nt = 6-12 lattices [66]. In the left panel of Fig. 1 1, trace anomaly 
is compared with those of p4 and asqtad quarks at ni„rf/mv = 0.05 on Nt = 8 lattices. The stout 
values are corrected by a tree-level improvement factor. We note that the stout data indicates 
slightly smaller Tc and much lower peak height of about 4 for (e — 3p) /T^. A generalized integral 
method was adopted to calculate the pressure in a multi-parameter system. They also attempted to 
calculate the EOS up to T = 1000 IVIeV. LCP up to large jS corresponding to this temperature was 
determined by a step-scaling method. Their results for the EOS are shown in the center and right 
panels of Fig. 11. 

4.2 EOS with 2+1 flavors of improved Wilson quarks 

About ten years ago, the CP-PACS Collaboration performed a systematic calculation of the 
EOS in two-flavor QCD with clover-improved Wilson quarks coupled to RG-improed Iwasaki glue 
[47, 67]. Extension of the study to 2+1 flavor QCD was not straightforward due to the heavy 
computational demand for Wilson-type quarks. 

To overcome the problem, Umeda et al. (WHOT-QCD Collaboration) developed a new ap- 
proach: a fixed scale approach combined with the T-integration method [68]. In the conventional 
fixed-A'; approach, the temperature is varied by varying coupling parameters on a LCP. In the fixed 
scale approach, on the other hand, the temperature is varied through a variation of A', with the 
coupling parameters fixed. The conventional integration method [65] in which the pressure is ex- 
pressed as a line integral of the trace anomaly in the coupling parameter space is inapplicable in the 
fixed-scale approach because the simulations are performed only at a point in the coupling param- 
eter space. Therefore, they developed an alternative integration method, the T-integration method 



P f 

[68] in which the pressure is given as an integration in T: ^ = / dT 



T5 
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The fixed-scale approach is complementary to the fixed Nt approach. Towards the high T limit, 
the fixed-scale approach suffers from lattice artifacts due to small Nt, while the fixed Nt approach 
is applicable for local quantities which are insensitive to the system volume. In the low T region, 
the fixed Nt approach suffers from coarseness of the lattice, while the fixed-scale approach keeps 
the lattice spacing a. With the fixed scale approach, a LCP is automatically followed. Since a 
common zero-temperature results can be used to renormalize thermodynamic observables at all 
temperatures, the cost for zero-temperature simulations can be largely reduced. We may even 
borrow zero-temperature configurations on the ILDG.. [See Ref. [69] for a different application of 
the fixed-scale approach.] 




Figure 12: EOS in 2+1 flavor QCD with clover-improved Wilson quarks adopting the fixed-scale approach 
[70] . (Left) Trace anomaly and its decomposition into contributions proportional to the beta function for j3 
and those for if's. The latter contribution is multiplied by —1. (Right) EOS obtained by the T-integration 
method. 

Adopting the fixed-scale approach, Umeda et al. attempted the first calculation of the EOS in 
2+\ flavor QCD with Wilson- type quarks [70]. They borrowed zero-temperature configurations 
generated by the CP-PACS-i-JLQCD Collaboration using the clover-improved Wilson quark action 
with non-perturbatively determined csw and the RG-improved Iwasaki gauge action [71]. Choosing 
the finest lattice with the lightest ni,,^, they carried out finite temperature simulations on A^f = 4-16 
lattices. The left panel of Fig. 12 is the result of the trace anomaly. They found a big cancellation 
between the contributions proportional to the beta function for jS and for K",,^ and Kg. The resulting 
low peak height is roughly consistent with those shown in Fig. 1 1 with highly improved staggered 
quarks, though the statistical error is large. Note that, the peak in the left panel of Fig. 12 is located 
at temperatures coiTcsponding to A'f 14. From the T-integration using a trapezoidal interpolation 
of the trace anomaly, the EOS shown in the right panel of Fig. 12 is obtained. The dominant part 
of the errors are due to a large cancellation in the zero-temperature subtraction of the gauge action 
in the low temperature region. Further statistics on low temperature lattices are needed to reduce 
the error bars. 

5. Further hot issues 

5.1 Chiral magnetic effect 

Two nuclei colliding with a small offset in heavy-ion collisions can create quite strong mag- 
netic field at the collision point [72]. When the gluonic topological charge fluctuates locally, the 
strong magnetic field may yield an electric charge asymmetry with respect to the reaction place 
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in the producted particles, as shown in the left panel of Fig. 13. This apparent violation of par- 
ity is called the "chiral magnetic effect" [73]. The STAR Collaboration reported an experimental 
evidence in favor of the predicted charge asymmetry last year [74]. 

Several lattice groups have studied the influence of external magnetic field on the quark matter. 
Buividovich and Kalaydzhan presented a study of quenched SU(2) QCD with valence overlap 
quark under a strong background magnetic field [75]. They observed an enhancement of electric 
current in the direction of the magnetic field, in accordance with the charge asymmetry. They also 
found that a strong external magnetic field induces nonzero electric conductivity in the direction 
of the external magnetic field in the low temperature phase. At the conference, D'Elia argued that, 
under a strong magnetic field, the transition temperature increases and the transition becomes 
sharper, as shown in the center panel of Fig. 13, based on a simulation of two-flavor QCD with 
unimproved staggered quarks at A^, = 4 [76]. 




Figure 13: (Left) Strong magnetic field and charge asymmetry created at the collision point of heavy ion 
collisions [74] (Center) Plaquette histogram as a function of the external magnetic field b in two-flavor QCD 
with staggered quarks on a 16^ x 4 lattice [76] (Right) Charmonium spectral function for the vector channel 
at zero temperature obtained on a quenched anisotropic 20^ x 160 lattice [82]. Blue symbols are the results 
of the variational method using different number of the basis functions. The red "+" are the results of MEM. 
Vertical dashed lines indicate the ranges of vector charmonia masses in nature. 

5.2 Charmonium spectral functions 

Nonaka studied charmonium spectral functions using the maximal entropy method (MEM) 
[77] on quenched anisotropic lattices with large spatial volume with and without momentum, and 
discussed the stability of Tj^ and above the transition temperature [78]. Allton calculated 
non-zero momentum spectral functions by MEM, using twisted boundary conditions on quenched 
isotropic lattices [79]. With staggered valence quarks, even time slice and odd time slice spectral 
functions are evaluated separately and combined to obtain the physical spectral function. From the 
longitudinal vector spectral function, they extracted the diffusivity. 

Dilepton production may be a key process to explore the formation of QGP. To calculate the 
dilepton rate and electrical conductivity without resorting to MEM, Francis and Karsch studied 
the vector meson correlation functions on large and fine quenched lattices in the high temperature 
phase using clover valence quarks [80]. They performed continuum extrapolation of them and 
calculated dilepton rate and electrical conductivity. Ding also presented their study of charmonium 
spectral functions by MEM, using clover valence quarks on quenched isotropic lattice [81]. To 
calculate a difference of spectral functions in which large energy parts are canceled and the low 
energy part we want to know is enhanced, they applied an "extended MEM" which is applicable 
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for not positive-definite spectral functions. Their result suggests that is melt and T]^ starts to 
dissolve at 1.467^. 

Ohno proposed a new method to calculate meson spectral functions adopting a variational 
method [82]. On finite volume lattices, the spectral functions should consists of discrete spectra 
only. Therefore, a precise meaning of the continuous spectral functions obtained by MEM is not 
clear. The variational method provides us with a method to directly evaluate the height and the 
location of low lying discrete peaks. In the right panel of Fig. 13, the heights and the locations of the 
spectral function by the variational method are compared with those obtained by the conventional 
MEM for the ground and first exited states in the vector channel on a quenched anisotropic lattice 
at r = 0. For MEM, the area of the peak is identified with the height of discrete spectra. The 
results for the ground state by the variational method fully reproduce the location and the area of 
the first peak obtained by the conventional MEM. For the first exited state, the variational method 
converges to a different point with increasing the number of the basis operators, which is closer to 
the experimental mass. Their preliminary result at T > indicates no cleai^ signs of dissociation 
for and rjc up to \ATc. 

5.3 Nature of hot quark matter 

Hydrodynamic evolution of quark matter plays an important role in the phenomenological 
analyses of heavy ion collisions. New attempts to calculate transport coefficients were presented 
by Kohno [83] and Maezawa [84] at the conference. 

Banajee studied scalar mesons around the transition temperature with two flavors of unim- 
proved staggered quarks on A', = 6 lattices with various spatial volumes [85]. They found that 
scalar meson, which is known to be unstable at T = 0, looks stable at T = 0.947^. Kitazawa stud- 
ied spectral properties of quai^ks in the Landau gauge in the deconfined phase on a quenched lattice 
with large spatial volume [86]. His result suggests a minimum at finite momentum in the plasmino 
dispersion relation. 

Three years ago, Garcfa-Garcfa and Osborn argued that the chiral transition of QCD may be 
viewed as a metal-insulator transition driven by the Anderson localization, based on a study low- 
lying Dirac eigenstates in quenched and 2-1-1 flavor QCD with 1-loop Symanzik gage action and 
and asqtad quarks [87]. A succeeding study of two flavors of unimproved staggered quarks by the 
Indo-French Collaboration concluded that the localization of low-lying eigenstates seems to be a 
finite volume artifact [88]. At the conference, Kovacs presented a new study of quenched SU(2) 
QCD and argued that the results of the Dirac spectrum are in favor of the Anderson localization 
picture [89]. More works are needed to clarify the status. 

5.4 QCD-like theories 

At the conference, Datta presented their study of the EOS in S\J{Nc) pure gauge theories for 
Nc = 3-6 [90]. They concluded that the speculated "strongly coupled conformal properties" (the 
diagonal dashed line in the plot) ai^e not seen. Similar result was obtained by Panero for Nc = 3-8. 
The A^c-dependence turned out to be mild, and, by adjusting the parameters, good agreement with 
a prediction of a holographic QCD model expected in the large Nc limit was discusses [91]. Feo 
presented a study of the EOS in 2-1-1 dimensions for Nc = 2-6 and argued that low-lying glueballs 
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do not explain the trace anomaly data in the low temperature phase [92]. For QCD-like theories 
with various Nc, many studies in different focusing have been put on the archive last year. An 
incomplete list includes [93, 94, 95]. 

6. Conclusions 

There were quite a few important advances in finite temperature QCD on the lattice in the past 
year. Big progress was made, in particular, with improved staggered quarks: The fruitful conflict 
about the value of the QCD transition temperature was (almost) settled. The 0(2) chiral scaling 
was finally observed by reducing the light quark mass down to a quite low value. The EOS and 
thermodynamic quantities are being measured in 2+1 flavor QCD around the physical point. 

From the intensive efforts to clarify the origin of the discrepancy in the transition temperature, 
we have learned the importance of the effect of taste violation in the simulations with staggered- 
type quarks. By correctly handling it — e.g. by adopting a highly improved action to sufficiently 
suppress the taste violation and by consulting the RMS hadron masses to judge the effective value of 
quark masses and to estimate the magnitude of taste violation effects — simulations with staggered 
-type quarks may achieve a quite high precision. 

On the other hand, the woiTy about the fundamental formulation of the staggered-type quarks 
remains. Therefore, it is indispensable to cross check the results with other theoretically sound 
lattice quarks. We have seen that steady advances are being made with Wilson-type quarks. Cal- 
culation of the EOS in 2-1-1 flavor QCD has started with improved Wilson quarks. Furthermore, 
simulations with the domain-wall quark action are now entering a level of quantitative investiga- 
tions. Lattice chiral quarks will open many applications beyond QCD. To compare with the results 
from the staggered-type quarks, however, further reduction of light quark masses towards the phys- 
ical point is necessary. 

I thank the organizers of Lattice 2010 for stimulating conference and their supports. I am 
grateful to D. Banerjee, A. Bazavov, P. Bicudo, S. Borsanyi, F. Bruckmann, M. Cheng, M. Chern- 
odub, S. Datta, M. D'Elia, C. DeTar, S. Ejiri, A. Feo, Z. Fodor, C. Gattiinger, F Karsch, T. Kovacs, 
A. Maas, H. Ohno, M. Panero, O. Philipsen, H. Saito, K. Szabo, G. Schierholz, W. Soldner and T. 
Umeda for sending me their new results and comments. I also thank the participants of Lattice 2010 
for useful discussions. This work is in part supported by Grants-in-Aid of the Japanese Ministry of 
Education, Culture, Sports, Science and Technology (No. 21340049). 
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